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Abstract 

Effects of large value assigned to the maximal car velocity on the 
fundamental diagrams in the Nagel-Schreckenberg model are studied by 
extended simulations. The function relating the flow in the congested 
traffic phase with the car density and deceleration probability is found 
numerically. Properties of the region of critical changes, so-called jamming 
transition parameters, are described in details. The basic model, modified 
by the assumption that for each car an individual velocity limit is assigned, 
is investigated in the aim to find the best supplementary rule allowing 
the jammed traffic to move with velocity larger than the slowest driving 
vehicle. 

PACS: 89.40. +k, 05.60.-k, 05.65.+b 45.70.Vn 

1 Introduction 

In recent years vehicular traffic problems have attracted much attention and 
a number of cellular automata models describing the traffic flow have been 
proposed in order to consider the dynamical aspects of the traffic system 012] ■ 
Presently, there are two basic cellular automata models that describe single 
lane traffic flow: the Nagel-Schreckenberg (NaSch) model |3] and the Fukui- 
Ishibashi model A, . Besides the cellular automata models, which are discrete in 
space and time, several other approaches to the traffic flow have been proposed 
recently. Among them there are space-continuous models in discrete time such 
as model of Krauss et al. or models continuous in space and time, e.g., 

the macroscopic (fluid-dynamical) models [7]. 



The cellular automata traffic modeling enables easy examination of micro- 
scopic and macroscopic aspects of the traffic flow and therefore it is considered 
as an important tool in traffic engineering • Though the models are fairly 
simple, it has been shown that they are able to reproduce real life traffic phe- 
nomena such as spontaneous formation of jams [8], explain the impact of global 
traffic light control strategies [§| on conditions for car accidents [TUUTT] . 

The basic NaSch model parameters are: a maximal allowed vehicle velocity 
"Umax, a vehicle density p, and a random deceleration probability p. The mean 
velocity v of all vehicles on a lane is the observable measured. At low density, 
p « 1, the vehicle flow is characterized by a linear dependence on the vehicle 
density. At high densities the vehicle flow decreases with increasing density and 
vanishes for the density p = 1. The vehicles are said to be in the congested 
state — so-called stop-and-go waves dominate in the dynamics of the system. 
Thus, with increasing car density the traffic state changes from the free flow to 
the congested traffic. 

The limit v rnax = oo has been introduced by Sasvari and Kertesz |12| . They 
have found that the fundamental diagram has a form quite different from that 
for a finite v max . The flow does not vanish in the limit p — > since even a single 
car produces a finite value of the flow. Due to the hindrance effect of other 
cars the flow is a monotonically decreasing function of the density. Hence the 
jamming transition disappears. 

In the following paper we investigate vanishing of the jamming transition. In 
particular, we study the influence of the general car velocity limit on the traffic 
features in order to identify the model parameters for which this transition 
disappears. 

Within our model several specific real life observations have been attempted 
to reproduced. Let us mention the German highway rules different rules 
for overtaking, see for the such rules review, or even the individual driver 
reaction to traffic ^3] . In this context we decided to consider a very Polish rule 
— each driver has its own velocity limit. One can give many reasons supporting 
our assumption: high diversity of vehicles, namely, from sport cars to scrap cars, 
or the common driver habit of breaking the speed limit. However, if at random 
we assign to each driver his/her own maximal driving speed then the stationary 
state is determined by the slowest vehicles. Therefore, again following hints 
from the Polish roads, the slowest driver is forced to change his/her behavior. 
By simulation, we observe the self-organization of the traffic when the maximal 
velocity of the slowest moving vehicle is changed. 

2 Standard NaSch model 

The probabilistic cellular automata model of a traffic represents a lane as one- 
dimensional lattice of L cells. Each cell is either occupied by one of N cars or is 
empty. A car can move with the velocity determined by integer values bounded 
by a speed limit v max - At time t a car is identified by a cell number ay which 
the car occupies and the velocity vf* with which the car go. The number of not 
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occupied cells before the car is denoted by gf' — xfh — x^p — 1 and is usually 
called the gap. The cars move along the lane following rules related with the 
driving habits: 

A Since every driver tends to drive with the maximal allowed velocity the 
car velocity is incremented in every time step (acceleration) 

V f+V =minK (t) + l,fw). (1) 

B If there is a risk of collision with the preceding car, the driver decreases 
the velocity to avoid the crash (braking) 

v\ =vom{v\\g\ J ). (2) 

C In order to include the unpredictable factors which influence reactions of 
the driver the car velocity is decreased with a small probability p (ran- 
domization) 

if random < p then vf +1 ^ — max(nf' — 1,0). (3) 
Finally, the position of the car is updated (movement) 
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Figure 1: The fundamental diagram of the NaSch model for v ma x = 5 and 
different random deceleration probabilities p. 

At each discrete time step t — > t + 1 the positions and velocities of all cars are 
adjusted in parallel, first, the velocities are obtained by the rules llll.'il) . then, the 
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positions of the cars are determined according to |@J. The periodic boundary 
condition is applied what effects in that the number of cars is conserved. 

Stationary states of cellular automata with the number conserving rule have 
in most cases a very simple structure which is reached after L - the lattice size, 
time steps ^H]- If J> = 0, since the evolution is deterministic, then the stationary 
state in the NaSch model at low car density is perfectly tiled by the two following 
patterns 

(e) and (e, e, . . . , e, v max ) (5) 

^ V 

where e means an empty cell and v max denotes a cell occupied by a car with 
the actual velocity v max . It is easy to check that this perfect tiling is possible 
only if the car density p satisfies the relation: 

P < l —T = Pent (6) 

If the car density p > p cr it then all intermediate to JSJ lattice patterns appear: 
(0), (e, 1), (e, e, 2), . . . , ( e,e . , e , v max - 1) (7) 

If the random deceleration parameter p > then all patterns from the list (J7J) can 
also be found for p < p cr u. The mean velocity is < v >— v max — P- Therefore, 
the parameter p can be considered as the symmetry breaking field and for the 
order parameter related to this field m — v max — < v >, the thermodynamics 
characteristics (critical exponents, scaling relations) of the second order phase 
transition can be estimated The similar results can be also found by the 
mean field approach |16| . 

In order to use the empirical results for a theoretical analysis it is convenient 
to use the mean flow J =< vp > for a given car density. The relation between 
the vehicle density p and the vehicle flow J is called the fundamental diagram. 
The flow depends on the car density and deceleration probability J(p,p). The 
famous fundamental diagram for v max = 5 is shown in Fig. ^ In this case 
Pcrit = 1/6 and J(p cr it,P — 0) is the maximal flow. When p is increasing then 
the densities which allows the car flow to be maximal decreases. In the following 
we will call the density corresponding to the maximal traffic flow as the critical 
car density for a given p and denote p cr it{v)- 

3 Simulation procedure 

We simulate a one lane traffic using the model described in the previous section. 
The one-dimensional lattice of the length L = 10 000 sites with the periodic 
boundary conditions forms the lane. If N denotes the number of cars on the 
lane then the car density p is defined as p = N/L. Our simulations start 
at random initial conditions; random localization of cars on a lane and random 
initial velocities. For each initial configuration, we update the individual vehicle 
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velocity and position in accordance with the update rules. The basic simulations 
are performed for the car density p e [0.01,0.5] with the step Ap = 0.01 and 
deceleration probability p 6 [0, 1] with the step Ap = 0.05. The results are 
obtained by averaging over 10000 time steps after the first 10000 ones left for 
the system stabilization. The procedure is then repeated for a number (100) 
of different realizations. The averages over different realizations provided the 
mean values of physical quantities which are presented in subsequent figures. 



4 Results 

4.1 Large velocity limit 

The limit v max = oo has been introduced in ^2] as equal to the length L of the 
simulated finite system. Here we propose to investigate large velocities, however 
still finite to observe results of the velocity limit to properties of stationary traffic 
on a lane. In the series of contour plots in Fig. |2]thc fundamental diagrams are 
presented when the velocity v ma x is increased from v max = 1 to v max — 500. 

If the velocity limit is low, e.g., w max < 10, then the fundamental diagram 
is a two part linear function, similar to that shown in Fig. ^ The first part 
corresponds to the linear increase of the car flow in case of free flow regime. 
The second part is the linear decrease that is characteristic for the congested 
traffic. There is a narrow intersection interval of p where the transition from 
one phase to the other takes place. When i> max grows then the transition point 
moves to lower densities. The congested traffic appears at all car densities with 
the exclusion of densities corresponding to almost empty lane, p << 0.1. 

It is noticeable that diagrams are almost indistinguishable when states of 
the congested traffic are considered for w max > 50. In Fig. |3|we present the 
numerical estimates for the linear decrease of the car flow J(p, p) in the congested 
regime for p > 0.2. The Pearson correlation coefficients for these estimates are 
r 2 > 0.99. Therefore, the global dependence of j on p and p can be found and 
our estimates lead to the following approximating formula: 

1 - 0.9p 1 - 0-8p 

J(M) — "W (8) 

It is easy to notice that in deterministic traffic limit, p = 0, we obtain the well 
known rigorous solution J(p, 0) = 1 — p. 

The properties which are dependent on the velocity limit occur at lower 
densities. Therefore, we focus on the free flow and especially on localization of 
the maximum of the traffic flow. Figs 0] and [5] are to show different properties of 
the free flow and the transition. Figure 01 presents all data obtained by us while 
in Figs EJa)-(d) the extracted data is shown in a way to visualize transition 
features. 

The traffic flow increases rapidly when both the car density and veloc- 
ity limit Vmax grow. The critical car density p cr it{p) and the maximal flow 
Jmax{pcrit(p),p) (see Figs G^c) (d) ) are related to the end of the linear increase 
of the flow in the free flow regime. 
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Figure 2: Traffic flow J for different global velocity limit v max = 
f, 5, 10, 50, 100, 500 - the contour plots. Notice, changes in the free flow re- 
gion and independence of the flow from the velocity limit in case of large v ma x 
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Figure 3: Simulation and numerical estimates for the expected stationary solu- 
tion of the traffic flow J = A{p) + B(p)p, independently of v max . 

Comparing the flow diagrams around the critical points presented in Fig. ^ 
to corresponding ones shown in Figs|S{a)(b)) one should notice that the value of 
the critical flow is significantly larger than it would be arisen from the congested 
state regime. One can suspect that we observe the occurrence of the free flow 
traffic in the congested state regime. Hence, possible that a kind of a regime 
with metastable states emerges for some interval of densities. 

Measurements on real traffic reveals that the traffic flow can exhibit metasta- 
bility and hysteresis effects but in the NaSch model in its simplest form these 
properties are not observed JTj. However, the metastable states and hystere- 
sis are observed in the so-called Velocity-Dependent-Randomization model |18| . 
Therefore it could be interesting to check the hypothesis whether starting the 
system with all cars stopped in a traffic jam one obtains a stationary state dif- 
ferent from a state reachable from a random initial condition. We performed the 
simulations for v max = 100 and it appeared that the stationary flows were only 
slightly lower than those presented in Fig. Ufa). The largest differences were 
observed in the interval 0.3 < p < 0.5 and they were less than 0.02. Thus, we 
cannot claim that in the NaSch model when the velocity limit is large a region 
with a hysteresis memory emerges. 

As it is expected, see Sec. 2 formula ©, when v max increases the pcritip) 
moves toward 0-density. Noticeable, that p influence to p cr it is weakened, see 
Efc) . For v max = 1000, p cr it = 0.001 for almost all p. If p > p cr it{p) then, as it 
is expected, the monotonic and independent of v max decrease of maximal flow 
is observed, see Fig. EJb) . However, the decay curve consists of the two easily 
recognizable parts: a nonlinear decrease (which can be approximated by some 
quadratic function with high accuracy) and linear decrease. The maximum of 
the flow J m ax also becomes p independent when the speed limit is sufficiently 
large, see Eld). 
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Figure 4: Traffic flow J for large global velocity limits v max — 50, 100, 500, 1000 
and low car density - the contour plots. The simulation step is Ap = 0.0005. 
Notice high peaks in traffic flows appearing in regions of a transition from free 
flow traffic to congested traffic. 

4.2 Individual velocity limit 

Let us assign the maximal velocity to each driver individually : v max (i),i = 
1, . . . , N. These velocities are chosen at random from the interval 1,2, ... , v max . 
Due to the hindrance effect of other cars J(p) is determined by the slowest 
moving drivers. Hence the obtained fundamental diagram is almost the same as 
the fundamental one obtained in basic NaSch model with v max = 1. It means 
that after stabilization steps all cars are moving in clusters which are lead by 
cars with the speed limit equal to 1. Therefore, we search for the supplementary 
rules that cause a faster traffic movement. 

let us start with the following rule: at each time step one of the drivers — 
that one who drives with the smallest actual velocity, is forced to revise his/her 
maximal velocity. Hence, before the acceleration substep (1), one car with the 
minimal actual velocity is found and then the new maximal velocity is assigned 
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Figure 5: Extracted properties of the transition from the free flow traffic to 
the congested traffic in case of large maximal velocities, (a) The fundamental 
diagram of the NaSch model for the same speed limit v max = 1 00 and different 
deceleration probability, (b) The fundamental diagram of the NaSch model 
for different speed limits v max and the same deceleration probability = 0.20. 
(c) The critical car densities obtained for different v max (log-plot), (d) The 
dependence of the maximal flow on deceleration probability, for different v max . 



to it at random. Let us denote this rule as (f , 0) while the basic NaSch model is 
denoted (0, 0). The resulting flows plotted for different p and p are shown in Fig. 
|SJ Their similarity to the shape of diagram with v max — 1 is still noticeable. 

To enhance the effect of the (1,0) rule, the following modifications are pro- 
posed (on the left, the notation of the new rules is placed): 

(2, 0): the new velocity limit assigned to one of the slowest car must be greater 
than the car previous limit 

(0, 1): if a gap between two subsequent cars is equal to 0, then the speed limit 
of the hindered car is increased by 1. The new speed limit cannot excess 
the overall limit v max . 

The combination of the above notations means that the corresponding combi- 
nation of the rules is applied simultaneously. For example, by (1, 1) we mean 
that one car with the minimal velocity is assigned a new velocity limit and all 
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Figure 6: Traffic flow when each vehicle has its own maximal velocity v max (i) 
chosen at random from 1, .... 10 (left panel) and from 1, . .., 100 (right panel). 
The slowest driving vehicle has the speed limit exchanged with a new one. 

cars which are separated from the following car by the zero gap have their speed 
limit increased by 1. 

Starting from random initial states with individual velocity limits taken at 
random from the interval 1,2, ... , 10, we observe the development of station- 
ary states. In Fig. the averages over a hundred realizations are presented. 
The time evolution of the mean velocity and changes in the average maximal 
velocity are shown for the two low car densities: p = 0.01 and p = 0.1, and 
low deceleration p = 0.05. The time axis is logarithmic to better visualize dif- 
ferent features of subsequent time intervals. In Figs. ©© we show the time 
development of small gaps = 0, 1, 2, 3, i.e. gaps responsible for the presence 
the following patterns (0), (e, 1), (e, e, 2) and (e,e,e,3), respectively, to have a 
deeper understanding of reasons for the global traffic properties. 

When the car density is low, p — 0.001, Fig[3a), then in 10 first steps all 
cars accelerate to move with their individual maximal velocities. The mean 
velocity < v > is only slightly lower than < v ma x >■ Only few small gaps are 
observed. Hence cars are separated suitable far away from each other on a lane. 
But then cars with low maximal velocity limits cause hindrance to a free move 
of faster driving cars. 

In (0, 0) model, i.e. in the NaSch model, the number of stopped cars and cars 
moving with < v >= 1 systematically grows while other velocities decay. If the 
slowest running car obtains a new maximal velocity at random, i.e., (1,0) rule 
is applied, then the increase of stopped cars and cars moving with < v >= 1 is 
slowed down what results in the average velocity higher than 2. However, after 
many time steps, both rules lead to the similar stationary state: more than a 
half cars drive with velocity equal to 1, and there is a large number of cars 
stopped. 

Surprisingly, that the mean maximal velocity < v max > settles at the value 
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Figure 7: Time development of the mean velocity < v^) > and mean maximal 
velocity < Vmax >) in simulations when each vehicle has its own maximal ve- 
locity v m ax(i) chosen at random from 1, . . . , 10. The labels at curves indicate 
at the supplementary rules applied to the basic NaSch model, (color plots on 
line) 



lower than the value obtained in case the (1,0) rule is applied. It happens 
because in our update configuration algorithm the slowest car is chosen as the 
most left of all cars moving with the slowest velocity. Hence, this car is often 
a last or close to the last car in a large traffic cluster. Usually, the tail cars 
are cars with the high speed limit. Therefore, a random change means here the 
exchange of the high speed limit to some another one, often lower. So that after 
many time steps we observe the decrease of the mean speed limit. 

If the slowest actually driving car obtains a greater velocity limit, i.e. the 
rule (2, 0) is supplemented, then < v max > must be a non decreasing function 
and, indeed, this function is growing fast in the first 5000 time steps. Due to 
the low car density, the number of slowly cars is small and the mean velocity 
is high. After many steps almost 1/3 of cars are separated by three empty 
cells: (e, e, e, 3), what allows the global traffic to move at the mean velocity 
< v >— 2.6 - the highest mean velocity observed in a stationary traffic in our 
experiments. However, if the density of cars is larger, see Figs db) and then 
the differences in traffic properties between the three rules considered become 
negligible. The slowly moving cars drive slowly because of hindrance effect, not 
due to own speed limits. 
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Let the NaSch model with random velocity limits be supplemented by (0, 1) 
rule. Then, it appears that the small number of stopped cars is permanent. The 
system stabilizes in a state where more than a half cars are separated by 2 cells. 
This effects in the greater mean velocity than in (0, 0) model, namely < v >= 2. 
Notice, that if between two subsequent cars the gap is 2 then two subsequent in 
time random brakes must happen to the first car to make the rule (0, 1) work, 
namely, to increase the first car maximal velocity. In case of p = 0.05 it is a 
rare event. The system stabilizes quickly, in less than 1000 steps, for both car 
densities simulated. 




H™ stop 

Figure 8: Time development of the small gaps — distances in empty cells be- 
tween subsequent vehicles. The labels at curves indicate at the supplementary 
rules applied to the NaSch model with individual speed limits. The vehicle den- 
sity is very low p = 0.01 and the deceleration probability is low p = 0.05. (color 
plots on line) 

If together with accelerating those cars which cause stopping others, the 
slowest running car obtains a new random velocity, i.e., the (1,0) rule is also 
supplemented, then a traffic preserves the property of having only low number 
of stopped cars together with low number of cars separated by 1 empty cell. 
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Figure 9: Time development of the small gaps — distances in empty cellsbctwccn 
subsequent vehicles. The labels at curves indicate at the supplementary rules 
applied to the NaSch model with individual speed limits. The vehicle density 
is rather low, namely p = 0.1 but higher than in the density in the xperiment 
presented in FiglHl The deceleration probability is p = 0.05. (color plots on 
line) 



In the time interval of 100, . . . , 1000 other properties of this traffic are similar 
to those observed when only (1, 0) rule is applied, though the mean velocity is 
higher. But in the stationary state the mean velocity becomes almost 2. In case 
of larger density of cars, the stabilization is reached faster and at approximately 
the same velocity. More than a half cars are separated by 2 two cells. The time 
development of systems governed by (0, 1) and (2, 1) is indistinguishable. 

5 Conclusions 

Despite its simplicity the NaSch model and its variants are able to reproduce 
many phenomena observed empirically. Although the same is true for other 
approaches mentioned here in the introduction, cellular automata model has 
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the big advantage of being ideally suited for large-scale computer simulations. 

However, the results presented in this article in the series of figures, due to 
high accuracy and minuteness of detail, are demanding a large computer time. 
But their clear structure ensures us in their reliability. 

In the following paper we basically investigated vanishing of the jamming 
transition when the velocity limit grows. Sasvari and Kertesz ^2] suggested 
that the behavior of the NaSch model is governed by the two fixed points p = 
and v max = oo. Schadschneider |16| indicated that one should rather consider 
p = and p = 1 as limiting fixed points in order to understand properties 
with finite velocity speed. Our results hint that for v max sufficiently large, the 
difference between the maximal possible flow, J max — 1 and J(p cr it is negligible 
small for all p £ [0, 1), i.e. for all p with the exception p = 1 where all cars 
are stopped from the beginning and the evolution cannot influence the initial 
car patterns. On the other hand, for any p ^ and p ^ the intermediated 
patterns of JZJl-list are present. The free flow regime vanishes. Hence, the limit 
Vmax — > oo excludes both deterministic evolutions: p = andp = 1. Properties 
of the flow in the congested traffic phase are independent of speed limit and 
only deceleration probability influences the value of the mean flow. 

By exchanging the lane velocity limit by the each driver velocity limit, we 
wanted to simulate the situation present on Polish roads. This situation results 
not only from the popular among drivers habit to break the speed limit but 
also from diversity of vehicles moving along the same lane: a horse wagon and a 
sport car. The simulations exemplify the dominant influence of the the slowest 
vehicle. However, the insistence of the car following the slowest vehicle , can 
improve, make higher, the velocity of the jammed traffic. 
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